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Abstract: The exact closed-form expressions for outage probability 
and bit error rate of spectrum sharing-based multi-hop decode- 
and-forward (DF) relay networks in non-identical Rayleigh fading 
channels are derived. We also provide the approximate closed-form 
expression for the system ergodic capacity. Utilizing these tractable 
analytical formulas, we can study the impact of key network pa- 
rameters on the performance of cognitive multi-hop relay networks 
under interference constraints. Using a linear network model, we 
derive an optimum relay position scheme by numerically solving 
an optimization problem of balancing average signal-to-noise ratio 
(SNR) of each hop. The numerical results show that the optimal 
scheme leads to SNR performance gains of more than 1 dB. All 
the analytical expressions are verified by Monte-Carlo simulations 
confirming the advantage of multihop DF relaying networks in cog- 
nitive environments. 

Index Terms: Ergodic Capacity, Rayleigh fading channels, Cog- 
nitive radio, Underlay relay networks, Amplify-and-Forward, 
Decode-and-Forward. 



I. Introduction 

Cognitive radio (CR) has drawn considerable attention 
in the academic and industrial communities in the past 
few years and has been considered as one of the main fea- 
ture for future wireless networks [TJH]. As an evolution of 
software-defined radio (SDR), CR with ability of learning 
from its surroundings and adapting its transmitting con- 
figurations allows secondary (unlicensed) users to oppor- 
tunistically transmits data in bands licensed to primary 
users [3H5]. As a result, it can alleviate the problem of 
spectrum congestion and thus allowing for a more efficient 
spectrum utilization [6]. 

Recently, CR has also been considered as the radio plat- 
form for relaying networks. Previous works on cognitive 
networks have assumed two types of cognitive operational 
modes including opportunistic spectrum access (OSA) and 
spectrum sharing (SS) [BJ[7]. In the OSA approach, CR 
users are allowed to transmit over the same frequency band 
licensed to primary users (PUs) only when the frequency 
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band is detected vacant [SrUOj. For the latter case, sec- 
ondary users (SUs) consisting of the source and relays may 
take advantage of a PU's frequency band by opportunis- 
tically transmitting with high power as long as the PU's 
transmission is strictly protected [H1[T2|. I n spite of the 
burden on PUs, the SS approach can improve spectral effi- 
ciency more aggressively than the OSA approach. 

Recently, the performance analysis for SS-based cogni- 
tive relay networks has gained great attention, (see, e.g. 
[13hT8] ). due to its adaptability to extend the coverage 
for wireless networks. This situation may arise in prac- 
tice when data is sent from a cognitive source to a given 
cognitive destination on a hop-by-hop basis via various in- 
termediate cognitive relay nodes. In addition to the advan- 
tage of extending the coverage without using large power at 
the transmitter, cognitive relay networks are able to reduce 
interference causing to PUs. 

In particular, Guo et. at in [13 derived the upper-bound 
for outage probability of underlay selective DF relay net- 
works operating within the constraint imposed on the peak 
power received at the primary receiver. In |14j . only based 
on the partial channel state information (CSI) of primary 
user link and partial CSI between the relays and the pri- 
mary user receiver, Li proposed two relay selection schemes 
for cognitive relay networks. The system performance in 
terms of outage probability was also provided over Rayleigh 
fading channels. For amplify-and-forward (AF) relaying, 
|15j and [16] studies the dual hop relaying networks with 
and without considering combining technique at the sec- 
ondary destination, respectively. Equipped multi antenna 
for secondary nodes, the work in [17] investigated the er- 
godic capacity for secondary underlay multi-input multi- 
output (MIMO) networks, where transmit antenna selec- 
tion (TAS) and maximal ratio combining (MRC) are used 
at the transmitter and receiver, respectively. In parallel, 
Dong Li |18j analyzed the effect of maximum-ratio combin- 
ing diversity on the performance of underlay single-input 
multi-output (SIMO) systems where the transmit power 
and the interference power constraint are taken into ac- 
counts. All of the theoretical performance analyses of cog- 
nitive relay networks mentioned above have just only been 
solved for the particular case of two consecutive hops except 
for the paper [19]. In [19], the authors studied the perfor- 
mance of cognitive underlay multihop decode-and-forward 
relaying networks over Rayleigh fading channels. However, 
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the results in this paper is limited while it only provided 
the closed-form expression of the system outage probability 
under interference constraints. 

Therefore, in this paper, we develop a performance anal- 
ysis framework for SS-based multi-hop CR networks. In 
particular, we derive exact closed-form expressions for out- 
age probability (OP) and bit error rate (BER) as well as the 
approximate closed-form expression for ergodic capacity for 
the considered system over independent but non-identical 
distributed (i.n.d.) Rayleigh fading channels. To gain in- 
sights, the asymptotic approximation for outage probabil- 
ity, bit error rate at high SNR regime, is also provided. We 
have shown that the system diversity order is always one 
regardless the number of hops and the coding gain increases 
according to the increase of the number of hops. For a pre- 
determined position of a primary receiver, the problem of 
relay position optimization is also considered and solved 
by using the numerical approach. Finally, simulation re- 
sults are provided to validate the analytical performance 
assessments. 

The remainder of this paper is organized as follows. In 
Section [TJ the system model for underlay cognitive multi- 
hop network is described. In Section lull the system per- 
formance metrics in terms of outage probability, bit error 
rate and ergodic capacity are derived for Rayleigh fading 
channels. In Sect. IIV1 we are concerned with the problem 
of relay position optimization. Numerical results are given 
in Section |V1 where the advantage of cognitive underlay 
multhop systems is investigated. Finally, conclusions are 
drawn in Section fVTl 

II. System Model 
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Fig. 1. System model of a multihop spectrum sharing communication 

system . 

We consider a multi-hop SS system with the coexistence 
of PUs, i.e., licensed users, and SU, i.e., unlicensed users, 
as shown in Fig. [T] Assume that PUs and SUs share the 
same narrow-band frequency with bandwidth B, which is 
licensed to PUs. In the primary network, the PU trans- 
mitter (PU-Tx) transmits its data towards the PU receiver 
(PU-Rx). In the secondary network, the SU-Tx (CRi) indi- 
rectly transmits the message towards the SU-Rx (CRk+i) 
with the help of K-l cognitive decode-and-forward (regen- 
erative) relays in between denoted by CR2, . . . , CRk ■ Fol- 
lowing the SS-based (underlay) approach, SUs are allowed 



to concurrently transmit with PU over the same licensed 
band while adhering to the interference constraint on the 
PU-Rx, i.e., any data transmission of SUs resulting in a 
higher interference level than an interference temperature 
at the PU-Rx is prohibited. To represent the maximum 
allowable interference power level at the PU-Rx, the in- 
terference temperature (Z p ), is used [T2TTPI] . Let hn,k 
and hi,k be the channel coefficients of the link from the 
fc-th SU-Tx to the next SU-Rx and to the PU-Rx, re- 
spectively. Under Rayleigh fading, \ho,k\ 2 and |/i/,fc| 2 are 
exponential distributed with their corresponding parame- 
ters X D ,k = E{\h D , k \ 2 } and \ Ijh = E{|/i J>fe | 2 }, where E{.} 
stands for the expectation operation. Perfect channel state 
information (CSI) of the CR& — > PU's Rx link is assumed 
at the CRaOJ We further assume that the additive white 
Gaussian noise (AWGN) associated with each hop is a zero- 
mean Gaussian random variable with variance A/"o . 

III. Performance Analysis 

Consider the fc-th hop of CR multihop networks, the in- 
stantaneous signal-to-noise ratio is written as 

2 



Ik 



Pk 



lD.k\ 



fl 



To ensure that the interference power at the PU-Rx is al- 
ways below the interference temperature, I pi the transmit 
power is upper bounded by 

Pk <Z P /\hi, k \ 2 . (2) 

Aiming to enhance the system performance, we adopt 

P k = l p /\h I , k \ 2 , (3) 

yielding 7fc = 

For Rayleigh fading channels, the probability density 
function (PDF) of 7z, fc = \h z , k \ 2 with Z e {D,I} is of 

the form f 1ZM (j) = , where ^ z ,k = A z ,fc. The 

PDF of the received SNR at hop k, jk, is derived as [351 P- 
187, eq. 6-60] 



A* (7) 



-4> 

{i+akY 



•\ T 

\ P 



(4) 



where otk — ^^T p /No- From ((U), the corresponding cu- 
mulative distribution function (CDF) is given by 



^(7) 



/ 7 „ {x)dx 



7 



(5) 
7 + ak 

Having the PDF and CDF of each 7^ in hands, we are 
now in a position to derive the performance metrics of the 
system including outage probability, bit error probability, 
and ergodic capacity. 

It can be realized by many mechanisms, e.g., direct feedback from 
primary receivers, indirect feedback from band manager 1201 or using 
CSI of the PU-CRfc link with channel reciprocity property 12 1 1 . 
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A. Outage Probability 

In an interference-limited multi-hop regenerative relay 
system, an outage event is declared whenever any of the 
instantaneous forwarded SNRs of the K hops falls below 
a given threshold, jth- In other words, the overall system 
outage is dominated by the weakest hop. Thus, the end- 
to-end (e2e) outage probability is written as [23] 



OP = Pr [min(7i, . . . , j K ) < j th ] . 



(6) 



Since 7fc with k = 1, . . . , K are assumed to be independent 
of each other, ([6]) is rewriten as 



K 



OP = l- [J[l-F 7fc ( m ) 



(7) 



fc=i 



Substituting (|S|) into ([7]), we have 

op = i - TT — ^— . 



(8) 



Theorem 1: At high SNR regime, the system OP can 
be approximated as 



K 



X 



I.k 



(9) 

Proof: We start the proof by using the fact that the 
cross-terms, F lk ('jth)F'n{jth) with k ^ I in (j?J), can be 
neglected compared to F lk (j t h) for values of interest. As a 
result, ([7]) can be approximated as 



K 



0P«^F 7fc ( 7th ). 



(10) 



fc=i 



5. 5ii£ Error Rate 

In this section, we study the most generalized scenario in 
a multihop network where all the single-hops in the route 
have the different statistical behavioiQ i.e., all the links are 
i.n.d. with different average channel power, 7 P ^7 9 . Taking 
into account the fact that a wrong bit transmission from 
node p to node q (q > p) is equivalent to an odd number of 
wrong single-hop bit transmission between both nodes and 
employing the recursive error relation, we have the exact 
e2e BER of the system as [2"5] 



K 

BER e 2 e = ^ ^ BERp 
P =i 



K 

n (i 

q=p+l 



2BER g ), 



(11) 



where BER p denotes the average BER for square M- 
ary quadrature amplitude (M-QAM) modulation (M = 
4 m , m = 1, 2, • • • ) in hop p and is given by 



oo 

BERp = J BER AWGN / 7p ( 7 )d 7 . 



(12) 



In the above equation, BERawgn is the instantaneous 
BER of hop p, namely [2H] 



BERawgn 



log 2 V M Vj 

E E ¥ n erfc (\^r) 

j = l n=0 



Mlog 2 VM 



(13) 



where Wj = (1 — 2~- 7 )vM— 1, u) n 



_ (2w+l) 31og 2 M 
— 2M-2 ' 



and 



(-1)' 



^= + |J)- Here, [-J and erfc(x) = 

OO 

are defined as the floor and complementary 

X 

error function, respectively. Substituting (|4]) and (|13|) into 
(|12p and swapping integration and summation order, we 
hav^| 



Plugging ([5]) into ([TP]) and making use x/(l +x) w a; for 
small £, we arrive at the desired result. This completes the 
proof. □ 

From (J9j) , it is worth noting that the outage probability 
at high SNRs is determined by the channel gain ratios be- 
tween the data and interference channels rather than the 
average channel powers. To gain further insights, we prove 
the following theorem. 

Theorem 2: The system diversity order and coding 

gain are G d = 1 and G c = ^y t h E X^T ) ■ 



Proof: Observing 



it is obvious to see that 
the system diversity gain is one and the coding gain is 

~,h E ^jj according to OP -> (G c I/AA )- Gd [24]. 

□ 

From Theorem [21 it is worth noting that similar to con- 
ventional multihop networks, the system diversity order is 
always one regardless of number of hops and the increase 
of hops results in an increase of the system coding gain. 



BER C 



log,VM Vj 

E 

j=l n=0 



E 



A/log 2 VM 



(14) 



where J v is defined as follows: 



Jp= I erfc (V^VT) 2 d l- 

(7 + otp) 



(15) 



Using integration by parts, it is shown that J v is of the 
form 



Up — 



7 + a v 



TT J 7 + OL % 
7=0 



dry. (16) 



Jx 



J, 



2 Other scenarios including independent and identical distributed 
(i.i.d.) case are a special case of the network under consideration. 

3 It should be noted that in this paper we only consider square 
A/QAM; however the employed approach could be easily extended for 
other modulation schemes such as MPSK, MPAM and rectangular 
A/QAM. 
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To determine Jp, we need to compute J\ and Ji. For J±, 
making use the l'Hopital rule, we have 



Jl 



erfc U/Uj/y) 



o. 



(17) 



For fa, introducing a change of variables, namely, u — y/y, 
we can rewrite J72 as 



J-2 



du 



= 2 



du - 



du 



(18) 



Together with the identities [13 eq. (3.23.3)] and [28] eq. 
(7.4.11)], it is shown that 

J p = 1 - ^U^e^ a " V^erfc (y^o^) . (19) 

Plugging (1T91) into ([14)) . we achieve the closed-form ex- 



pression for BERp as 



log 2 V MVj 

£ [l- v /w^e w " Q! P v /7rerfc(^/a^a^)] 
BER„ = J=1 "=° ^ ^ . (20) 



/ Aflog 2 VM 

Substituting (I2TJ1) into (fTTj) yields the average e2e BER. 



For i.i.d. fading channels, i.e., {ot p } p=1 — a, (jlip simpli- 
fies as (|21[) shown at the top of the next page. 

Theorem 3: At high SNR regime, the end-to-end BER 
of cognitive underlay multihop DF relaying networks oper- 
ating in Rayleigh fading channels is approximated as 



BER 



c2c 



K 

tS- Y] — , i.n.d. channels 

p=l 
Kg 
2ba ' 



(22) 



i.i.d. channels 



whom n - V^M-l j u _ 31og 2 M 

Proof: Observing (fTTT) and using the fact that at 

K 



high SNR regime the product term, fj (1-2BER,,), 

q=p+l 

approaches to one. We are able to approximate the end- 
to-end BER as 

K K 

= ^ BERp Yl (l-2BERg) 

p=l q=p+l 



BER e 



K 

^ ' BERp. 

p=i 



(23) 



By neglecting some of the higher order terms in (jT3j) , BER p 
is expressed as [26] eq. (18)] 



BER 



a erfc 



-d-y, 



(24) 



where a 



/M-i 



and b 



31og 2 M 
2(M-1) ' 



Employing the 



/Mlog 2 vm 

same steps as for (|2T)|) and then making use the infinite se- 
ries representation for the error function, i.e., eric(y/x) ~ 
(l — 35) for large x [29], we can have the desired re- 
sults as in (|2"2"|) . 

For i.i.d. case, it follows immediately from the result of 
the i.n.d. case by letting a = {ce p }^ =1 . n 

C. Ergodic Capacity 

Besides outage and bit error probability, the ergodic ca- 
pacity is another important performance measure, defined 
as the expected value of the instantaneous mutual informa- 
tion between the cognitive source and the cognitive desti- 
nation. The ergodic capacity C (in bits/second) per unit 
bandwidth can be expressed as 



C=- 



1 

K 



log 2 (l +7)/ 7e2c (7)d7, 



(25) 



where r y e 2 e is the equivalent instantaneous e2e SNR for the 
multi-hop CR network. To derive the ergodic capacity, we 
first need an expression for the PDF of 7 e 2e- However, with 
regenerative relaying, an exact closed-form expression for 
the PDF is not mathematically viable. For mathematics 
tractability, we use the approximation approach. Accord- 
ing to [30] , regardless of the modulation scheme used, r ) e 2 e 
can be tightly approximated as 



7e2e 



7e2e 



mm 7 fe . 
k=l,...,K 



(26) 



Having been widely adopted in the performance studies 
of DF relay networks (see, e.g., [HTJOU]), the advantage of 
this analytical approach is able to provide a mathematics 
tractable form for the CDF and PDF of the end-to-end 
SNR. Consequently, the PDF of 7 e 2 e is given by 



f%2e (7) = 



d~; 

E-M7) n (i-^M), (27) 

k— 1 n—l.n^k 



K 



where i^ e2e (7) = 1 II I 1 ~ F jk{l)]- Substituting (gj and 
© into ([27]), we get 



fc=i 



K 



f%2. (7) = 



(28) 



With the current form of ([2"B")) . it is very difficult to ob- 
tain the closed-form expression for the end-to-end capac- 
ity. To facilitate the analysis, we sort and renumber au 
in the ascending order as ai = • • • = a n = ft\ < 
■ ■ ■ < ctr 1 +r 2 +---+r N _ 1 +i = ■ ■ ■ = a ri +r 2 +-+r N = Pn and 



K with T]~ being a positive integer. Stated an- 



other way, /3i, • • • , Pn are distinct elements of a\, • • ■ , ctK- 
Using the partial-fraction expansion, (|28|) can be rewritten 
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BER„2„ — — 



1-1- 



Mlog 2 VM 



log 2 vM Vj 

E E 

j = l n=0 



K 



<jp n (l - ^w„ae"" Q \/7rerfc (^a)) 



(21) 



as 



A r„ 



/^(7) = nL«MEE 



Ar. 



«=1I=1 (7 + /3n) 



i+1 



(29) 



For the special case of /3 n = 1, (f36|) simplifies to 
1 f°° dj 



where A n j is the coefficient of the partial-fraction expan- 
sion determined as [331 134^ 



lbx2j ( 7 +l) 
1 



z+i ' 



Z 2 ln2' 



(37) 



A, 



Q{r n -l) 



For i.i.d. case, from (|27l) . we have 



Plugging d29j in ([25") . we get 



[(7+^r/^ 2c (7)] 



(31) 



/7^(7)=^[l-^ fc (7)] A "V 7fc (7) 



7=-/9« 



nil' A r„ 



(32) 



n=l (=1 

where Ii(fi n ) with / > 1 is the auxiliary function defined as 



log 2 (l+ 7 ) J 
l+i a r- 



la (fin + 7) 
Using integration by parts, we have 

log 2 (l + 7) 



(33) 



7=0 



I In 2 



o (l + 7)(7 + A.) 



(34) 



When / is an integer, after using partial fraction expansion, 
we have (|35| as shown at the top of the next page. Note 
that the integral in (|35|) is not converging due to the first 
term. To deal with the problem, by appropriate rearrange- 
ments and then performing the integrations, we obtain the 
final closed- form expression for Ii(fi n ) as (|36p. 

4 For convenience, the coefficients A n ; can be obtained more easily 
by solving the system of K equations, which is established by ran- 
domly choosing K distinct values of 7 but not equal to any n 1341 . 
Let us denote K chosen values of 7 as B u with u = 1, • • • , K, we can 
obtain a linear system of equations as 



EE 



r = E 



rijrj (7 + /?n) I+1 £l(7 + o*)lEU(7 + an) 



(30) 



A = C 1 D, where [.] T denotes the transpose operator; C 



a K X K matrix, whose entries are C u ,v 



1 



p-1 

v = q+ E r m ; D = [ Di ■■■ D u ■■■ D K ] T with D u 

Tfl = \ 

K 



with 



FL,— 1 I p 1 T E ? o ! \ an d u > 

(Bu+otn) ^-> 1 (B u +a k ) ' 



(7 + a) 



K+l ■ 



(38) 



Combining (j38|) and (|25|) . we have the e2e ergodic capacity 
for this case as 

00 

n „JC f lo S2(l + 7) 



(7 + ") 



A'+l 



a K I K (a). 



(39) 



Plugging ([55 ]) (or (|37 j) ) into {32) (or (J3S)), we have the 
closed-form, integral-free, expression for the system capac- 
ity. It is worth noting that our suggested method is precise 
and tractable with the determination of the appropriate 
parameters being done straightforwardly. Additionally, as 
C is given in a closed-form fashion, its evaluation is instan- 
taneous regardless of the number of hops and the value of 
the maximum interference temperature. 

It is of interest to compare the ergodic capacity of under- 
lay DF and AF multihop networks. The following theorem 
is provided to answer such the question. 

Theorem 4: For the same network and channel set- 
tings of underlay multihop networks, DF relaying provides 
slightly better ergodic capacity than its AF counterpart. 

Proof: Denote Cdf and Caf be the ergodic capacity 
for the underlay DF and AF multihop network, respec- 
tively. According to the min-cut max- flow theorem [35], 
namely the end-to-end system capacity cannot be larger 
than the capacity of each hop, the end-to-end ergodic ca- 
pacity for underlay DF multihop networks can be written 
as 



Cdf < min(ci,c 2 , . . . ,c K ), 



(40) 



where Cfc with k = 1, . . . , K is the Shannon capacity of hop 
fc, given by 



(41) 
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l In 2 



(/? n -l) J ( 7 + l) 



E 



{(fir, 



1)V 



a+i-k 



(35) 



IlWn) 



1 

/In 2 
1 

/In 2 



n -l) z /(7+T ~7 + /3„) d7 ~^ Wn-lfl 



log/3„ 



(7 + ft.) 



z+l-fc 



i-1 

E 



(36) 



Making use the Jensen's inequality, we easily see that 



1 

K' 



Cdf > -^E 7l ,... j7K {min [log 2 (l + 71), . . . , log 2 (l + 7^)]} . 

(42) 



Since the binary logarithm is strictly concave, Cdf can be 
rewritten as 



Cdf > ■^E 7 i,...,7k 



{log 2 [l + min(7i,...,7K)]}. (43) 



Based on the results reported in 



, i.e. 



1 



C AF < -h^ u ...nK {log 2 [1 + min( 7 i, . . . , lK ))} , (44) 



we can complete the proof after making a comparison be- 
tween (03J and flU). D 
Before moving on to the next section, here we would 
like stress that although DF performs slightly better er- 
godic capacity, it needs more complicated implementation 
as compared to AF [23] , 

IV. Relay Position Optimization 

In this section, we focus on the problem of relay position 
optimization. In particular, for given underlay DF multi- 
hop network parameters including the coordinates of the 
primary receiver, the secondary source, secondary destina- 
tion, and the number of hops, our problem is to find optimal 
positions for relays, which makes the system performance 
(in terms of outage probability or system error probabil- 
ity) minimize. We hrst provide the solution for the case of 
outage probability. 

For simplicity, we consider the network scenario illus- 
trated in Fig. [2j where all secondary nodes are ordered in 
the sequence and positioned along the straight line connect- 
ing the secondary source and the secondary destination. 
Such a model is mathematically tractable and well-adopted 
in the literature in studying multihop networks [37,38 . In 
addition, it is readily extended to the more generalized case 
of two dimension (2-D) networks. Interestingly enough, al- 
though the linear multihop network model is slightly sim- 
plified model of the real world, we can find it in practical, 
e.g. the communication between cars on a highway, or the 



01,(0,0) 




CR(1,0) 



Fig. 2. Cognitive underlay 3-hop DF relay network in a straight line. 



communication between road-side units placed along the 
road [MED]. 

We further assume that the overall distance between the 
source and the destination is normalized to one, i.e., 



ZD,1 



ID,K 



1. 



(45) 



where dn,k denotes the physical distance of data hop k. Un- 
der a predetermined position of the primary receiver and a 
fixed number of hops K , the problem of finding the opti- 
mal position of the relays that minimizes the system outage 
probability can be mathematically stated as follows: 



min ^ T, 4^ 

K 

d D.k = 1 

k=l 

d D ,k > 0, k = l,...,K 



(46) 



subject to 



Based on a single-slope distance-dependent path loss 
model for the average channel powers |41j . we can write 



A 



I.k 



A 



D,k 



di,k n _ 

dn k - V 



( ^£21 



I.k 



(47) 



where 77 > 2 denotes the path loss exponent. Recalling that 
r] typically has value of 2 in free-space environments and 
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up to 5 and 6 in shadowed areas and obstructed in-building 
scenarios, respectively [HI Table 4.2]. Combining (|46l) and 
(|4"7| and noting that jth ancl jf- are the given constraints, 
the optimization problem in (1461) can be simplified as 



K 

E 



subject tu { £/ D ' fc 

dc,fc > 0, k 



1. 



(48) 



A 



Theorem 5: For a given coordinate of primary re- 
ceiver (xp,yp), the secondary network outage probability 
achieves its minimum at 



OP r 



7th 
x 



K 



* n 



D,fc 



1 



(49) 



where <i^, k with fc = 1, ...,K are the optimal distance 
of hop fc, being the roots of the nonlinear system of A 
equations as follows: 



IDA 



(x P - (Ida) - (x P 2 



_ 2 



1 



VP 



K-l 

X P - d D.k 
k=l 



And d} k with fc = 2, 
relationship 

"D,l 



x P " 



)fe) 

«p 2 ) 



?7P 2 = 



yp 2 = o 



(50) 



. . , A are determined by using the 

d* 



l D,k 



rl* 

a D,K 



d *i,i 



(51) 



d*i t k d* IK 
Proof: To solve the above optimization problem, we 
can use Cauchy theorerrj^]. In particular, from (|48[) . we can 
have 



K 

E 

k=l 



dD,k 
dlk 



>K1 



lD,k 



(52) 



Equality holds if and only if 



d 



D,l 



d 



D,k 



d iA 



*i,k 



l D,K 



(53) 



Denoting d* D k and dj k as the optimal values making the 
equality occur, the minimized OP is 



7Jo Vfe=i x ' k . 



(54) 



5 In some mathematics books, this theorem also is named as the 
inequality of arithmetic and geometric means. 



To determine d* D k and d*j k , observing Fig. [2] and making 
use the Pythagorean theorem, we have 



di,k = yp z 



xp 



at-i 

^k=\ 



lD,k 



(55) 



for k = 2, . . . , K. It is obvious for the case of = 1 that 
diA = y/xp 2 + yp 2 . Using (|53"|) . (1551) is rewritten as follows: 

(56) 



xp 



fc-i 

E< 

p=i 



yp 



[xp 



yp 



iD.k 



IDA 



Combining (|45|) and (|56|) . a system of A equations for 
do,i, ■ ■ ■ , do,K is formulated as follows: 

dD,i H h rfzj.K -1 = 



^ + (x P -rf D , 1 r-(xp 2 + yP 2 )(^) =o 



yp 2 + (x P -tdD,k) -(V+2/p 2 )(^) =0 

n \ h — 1 / 

(57) 



With the current form of (|57j) . it seems impossible to obtain 
the closed-form expression for do,k- Consequently, in this 
case, the only possibility is to solve (|57l) numerically. Using 
Newton's method [43 , (|57j) can be determined by means of 
the recursion. The initial value for da,k can be selected as 
d*D fc(0) = Tc f° r au From ([53)) and ([57jl. we can complete 
the proof. □ 
We next present optimal relay positions, which minimize 
the average system bit error rate. Different from outage 
probability, which serves as a lower bound to the frame 
error rate for block fading environment and provides an in- 
sight into the theoretic-information performance limit, the 
average bit error rate shows the actual system performance 
for a desired target spectral efficiency, i.e. modulation level. 
As such, the following theorem is of importance in this re- 
gard. 

Theorem 6: For a predetermined coordinate of pri- 
mary receiver (xp,yp), linear DF multihop networks un- 
der interference constraints provide the best performance 
in terms of bit error probability if and only if 



"■da 
d i,i 



l D,k 



d I, 



U D,K 
d I,K 



(58) 



And the corresponding system bit error probability under 
optimal relay positions is 



A' 



BER e2e = — A ( Y[ 



l D,k 



(59) 

Proof: The proof is omitted here due to the similarity 
of the form between OP and BER at the high SNR regime. 
Then the Theorem [6] is easily inferred from Theorem [5] □ 
From Theorem [5] and [51 it is worthy to point out that the 
optimal relay positions (both in terms of outage probability 
and bit error rate) do not depend on the path loss exponent 
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V. Numerical Results and Discussion 

The purpose of this section is twofolds. We first provide 
numerical results to confirm the derived analytical expres- 
sions and then show the network performance advantage 
offered by relay position optimization. 

For illustrative purpose, we consider a linear multi-hop 
network in a 2-D plane, where all SUs are co-linearly lo- 
cated and the distance between the cognitive source and 
the cognitive destination is normalized to one. Further- 
more, the cognitive source and the cognitive destination are 
located at points with coordinates (0,0) and (1,0), respec- 
tively. Each cognitive relay node is equidistant from each 
other, i.e. dcR fc ,CR fc+ i =1/K. The average channel power 
for the transmission between node A and node B is modeled 
as Aa,b = g?a,b~'' where r\ denotes the path loss exponent 
with A e {CRi, . . . , CRk-i} and B e {PU, CR 2 , . . . , CR K }. 
In all examples, we locate the PU-Rx at coordinate (0.35, 
0.35) and set r\ — 4. 



10 



10 



10" 



10 



10 





16-QAM 




/(decreasing 5 -» 1 


Analysis (Exact) 




Analysis (Approximation) 




O Simulation 









-10 



10 



20 



30 



40 



/ p /W [dB] 



Fig. 4. Bit error probability versus average Xp/Mo- 




Fig. 5. Ergodic capacity versus average Xj,/A/"o- 



Fig. 3. Outage probability versus average X p /A/b- 



in Fig. [3] and Fig. [4J we respectively illustrate the outage 
probability and BER of the multi-hop cognitive networks 
as a function of interference temperature for different num- 
ber of hops. As can be observed from the two figures, the 
performance is enhanced as the number of cognitive relay 
hops K increases. It is important to note that the dimin- 
ishing gain returns as the number of hops increases. The 
numerically evaluated results demonstrate the correctness 
of the presented analysis. 

Figure [5] displays the capacity performance for cogni- 
tive multi-hop transmission by varying the number of hops, 
K = 1,2,. ..,5. It is worth noting that for interference- 
limited regime (low T p /J\[q), the system with large K offers 
improved performance. For a high interference tempera- 
ture level, multi-hop transmission with small hops is more 
favorable. It can be explained by using the fact that with 
the channel model and time-sharing schedule used, at low 
Xp/No transmission over shorter distance corresponds to 
increased effective SNRs while at high X p /J\fo increasing 



the number of hops is equivalent to reducing the effective 
transmission bandwidth of each hop. In addition, the nu- 
merical results show that the analytical results are in good 
agreement with the simulation results. 

Fig. [5] compares the capacity performance of multi-hop 
cognitive relay networks for different positions of the PU- 
Rx given the same number of hops. Observing the results 
in the figure, we can see that the system performance im- 
proves when the primary node is located farther away from 
the secondary relay transmitters, as expected. 

Up to this point, we have not studied the effect of the 
proposed relay position optimization. In doing so, we con- 
siderer three relay position profiles: randomization, equal- 
ization and optimization, denoted as Profile A, Profile B 
and Profile C, respectively. In Profile A, all secondary re- 
lays are chosen randomly from a uniform distribution. In 
Profile B, the distance between any two nodes is the same. 
And in Profile C, secondary relays are set using the rule, 
proposed in Sect.[lV] Table Q] demonstrates the results for 
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Fig. 6. Effect of PU-Rx's location on ergodic capacity. 



Fig. 8. Effect of path loss exponent on ergodic capacity. 



K = 2 and 4. 





Profile A 


Profile B 


Profile C 


K = 2 


d D>1 = 0.1767 
d D ] 2 = 0.8333 


d D .i = 0.5 
d D ',2 = 0.5 


d DA = 0.4192 
rf D , 2 = 0.5808 


K = 4 


d D ,i = 0.20 
d D ,2 = 0.28 
d Dy3 = 0.36 
d DA = 0.16 


d D ,i = 0.25 
du,2 = 0.25 
d D ,3 = 0.25 
d DA = 0.25 


d D ,i = 0.1915 
efo, 2 = 0.1900 
d D [ 3 = 0.2492 
d DA = 0.3693 



Table 1. Comparision of three relay position profiles 




Fig. 7. Comparison of three relay position profiles. 

In Fig. we can see that Profile C outperforms Profile 
B, which, in turns, outperforms Profile A. This observation 
also repeats in Fig.[8l where the effect of path loss exponent 
is investigated. From Fig. [8j we can see that the optimal 
number of hops is a complicated function of rj. Further- 
more, the advantage of profile C, i.e. the capacity gap, 
becomes bigger with higher value of rj. With small 77, the 
increase of hops results in the loss of the ergodic capacity. 
However, with large values of 77, there exists a value of K, 



that makes the system ergodic capacity maximize. It can 
be explained by using the fact that with small 77, the ben- 
efit of path loss gain is not enough to compensate the loss 
due to the use of multi orthogonal time slots for multihop 
communications. 

VI. Conclusion 

We have investigated the performance of cognitive re- 
generative multi-hop relay networks using the underlay ap- 
proach. We have derived the closed-form expressions for 
the outage probability, BER, and ergodic capacity over 
i.n.d. Rayleigh fading channels. High analysis for outage 
probability and bit error rate have also made to provide 
insights into the system behaviors. The numerical results 
show that under the interference constraints inflicted by 
the primary network, the multi-hop transmission still of- 
fers a considerable gain as compared to direct transmission 
and thus makes it an attractive proposition for cognitive 
networks. 
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